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基于 互补 约束 规划 模型 的 投影 Barzilai-Borwein 
梯度 算法 求解 绝对 值 方 程 
严 涛 
(南京 理工 大 学 理学 院 数学 系 , 南京 210014) 


摘 要 : 在 绝对 值 方程 4x- 几 =2 问题 有 解 情形 下 ， 给 出 了 求解 绝对 值 问题 的 一 新 方法 。 首 先 建立 了 一 等 价 求解 绝对 
值 问题 的 互补 约束 规划 模型 ， 进 而 利用 新 模型 中 的 非 负 约束 ， 给 出 了 求解 绝对 值 问题 投影 Barzilai-Borwein(BB) 梯 度 
算法 。 数 值 实验 结果 表明 了 方法 的 有 效 性 。 
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Projected Barzilai-Borwein gradient method for absolute value equations based on 
mathematical programs with complementarity constraints model 


Yan Tao 
(Dept. of Mathematics, College of Science, Nanjing University of Science & Technology, Nanjing 210094, China) 


Abstract: This paper presented a mathematical program with complementarity constraints (MPCC) model to solve the 
absolute value equation Ax-|x|=b . Based on the new model, the projected Barzilai-Borwein (BB) gradient method was 


applied to solve a penalty subproblem and the solution of the original problem was obtained. Preliminary numerical results 


illustrated the efficiency of the new method. 
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其中 4=(4-7-4- 门 。 


加 7 定理 1 如 果 x* 为 绝对 值 方程 问题 (1) 的 解 ， 则 
J 考虑 如 下 绝对 值 方程 问题 有 
QQ Ax-IN=b (区 本 (区 妆 ] 为 昌 题 (2 的 解 : 反 之 ,如 果 w 为 


其 中 :， AeR”,beR" ，|x| 按 x 分 量 的 绝对 值 获取 。 绝 对 值 问 

题 为 广义 绝对 值 问 题 4x+B 思 =2 的 特殊 情形 吕 , 同时 在 一 定 条 ”问题 2) 的 解 ， 定 义 忆 = 呆 - 场 。 (=42.…, 有 ) , 则 x 为 绝对 值 方 
件 下 ,问题 等 价 于 线性 互补 问题 C3。 绝对 值 方程 问题 已 吸引 ” 程 问题 (1) 的 解 。 

许多 学 者 对 其 理论 和 算法 进行 深入 的 研究 ， 且 已 有 许多 相关 证 明 将 w 直 接 代入 问题 (2) 则 得 到 第 一 部 分 结论 。 对 于 


人 。 成 果 出 现 。 文献 [4] 对 相关 理论 进行 了 详细 介绍 。 基于 非 光滑 。 | y 
明 第 部 分 ， 首 将 ww 分 解 类 7 ,WT ”其 
方程， 文献 [5.6] 提 出 了 相应 的 半 光 滑 牛 顿 法 ， 应 用 各 类 光滑 。 第 一 部 分 首先 将 ww 分解 为 ww") 下 


函数 将 原 问 题 转换 为 光滑 方程 , 文献 [7~9] 提 出 了 相应 的 光滑 WE=W = ， 则 问题 (2) 中 的 第 一 个 等 式 可 
化 方法 ; 文献 [10] 在 4 为 对 称 和 矩阵 情形 下 ， 将 绝对 值 方程 问 写成 401 -两 )-(W' + 更 )=b5 。 因此 ， 可 构造 x=W'- 殉 ， 即 
题 转换 为 了 无 约束 优化 问题 ， 进 而 提出 了 求解 原 问题 的 = 坟 一 Wn(i=42,…,n) 下 面具 需 证 |xzFW+ 鸯 。 由 问题 (2) 有 
Gauss-Seidel 算法 。 另 一 方面 ， 互 补 约束 数学 规划 是 一 类 约 =0,ww, =0， 则 

束 中 含有 互补 问题 的 优化 问题 。 它 是 均衡 约束 规划 的 一 个 特 
殊 子 类 ， 且 在 工程 设计 、 经 济 建 模 等 实际 问题 中 有 着 很 好 和 
应 用 00。 近 年 来 ， 对 于 互补 约束 规划 的 算法 已 有 很 好 的 研究 


Pe ee = NW + wi) = W + wi 
展 ， 如 罚 函 数 逼 近 法 [9 光滑 逼近 法 0 正则 方法 上。 受 上 


er 


WE 


述 结果 启发 ， 本 文 建立 了 求解 绝对 值 方程 问题 的 互补 约束 规 。 则 |xEW+WW ， 第 二 部 分 结论 得 证 。 

划 等 价 模型 ， 进 而 构造 了 罚 函 数 子 问题 ， 并 用 投影 BB 梯度 进一步 ， 本 文 将 问题 (2) 改 写 为 如 下 互补 约束 规划 问题 : 
方法 求解 子 问题 以 期 获得 原 绝对 值 方程 问题 的 解 。 ee 
、 mn 

1 互补 约 束 规 划 模型 st. G(w)>0,H(w)>0, (3) 
首先 给 出 一 个 与 绝对 值 方程 问题 等 价 的 约束 方程 系统 ， GOO HOW SY, 

即 寻 找 we R* 使 得 其 中 :G0w)=w,Hi(W)=wis(i=14,2,.…,n) 。 且 w 为 问题 (3) 的 最 优 


汉 解 当 且 仅 当 ”为 问题 (2) 的 解 。 因 此 ， 求 解 绝对 值 方程 问题 
(1) 可 转换 为 求解 互补 约束 规划 问题 (3)。 


I 


Aw =b, 解 解 


2) 


Ww >0, ww =0,1=b,2,.,n, 
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Po gd) 计算 so =ze-zo ， y?=Vf(ztD) -Vf(z0) 。 如 果 
2 算法 


《50,y0)<0， 令 a8? =10”; 否则 令 

算法 包含 两 个 迭代 过 程 ， 一 是 用 来 更 新 变量 wo 的 外 迭 cl =min 10, max ho™, C0) | 
代 ， 另 一 个 是 用 来 求解 每 步 中 子 问题 的 内 欠 代 。 为 实现 求解 《07 
法 ， 由 文献 [15], 首 先 构 造 如 下 罚 函 数 问题 


6) 令 1=1+1l， 转 步骤 b)。 
(4) 文献 [16] 中 定理 2.2 有 如 下 结 
人 定理 3 ”由 内 迭代 算法 产生 的 序列 {zo} 的 任何 聚 点 都 是 
其 中 :Pp 为 罚 参量 。 模型 (4) 具 有 互补 约束 规划 问题 (3) 中 的 。” 罚 函数 问题 (4) 的 稳定 点 。 
简单 非 负 约 束 ， 且 可 方便 应 用 投影 算法 ， 同 时 本 文 可 以 证 明 。 。 数值 结果 
通过 求解 罚 函 数 问题 (4) 可 获得 互补 约束 规划 问题 (3) 的 稳定 
点 


和 
min zl Aw—bll + pG(w)’ H(w) 


。 因 此 ， 算 法 构造 将 以 模型 (4) 为 基础 。 外 和 迭代 算法 中 的 停机 准则 设 定 为 14zx-|zl-il<10* 。 

下 面 给 出 外 迭代 算法 过 程 及 相应 收敛 结果 。 下 面 讨 论 内 迭代 算法 的 停机 准则 。 ”对 于 简化 的 罚 函数 
1) 外 和 迭代 算法 模型 (5)， 本 文 定义 f(w) 投影 梯度 为 
a) 设 置 初始 点 w® = ,ww 20，P,B>1， 令 k=1; (worowm) = (Vf (Ww)), w>0 
b) 令 =w we,i=42,…,n， 如 果 1Ax-|x|-bl=0， 停 止 计 ?lmin{0,(Vf (Ww)),} w=0. 

算 ; 则 w 为 问题 (5) 的 稳定 点 当 且 仅 当 WV?f(w)l=0 。 另 一 方面 由 文 
0) 令 p=Pi， 以 ww 为 初始 点 , 求解 问题 (4) 的 稳定 点 ww; ” 献 [16] 的 引 理 2.1 可 知 , 上 面条 件 等 价 于 1PCwW-V0w)-wll=0， 
中 ) 令 k=k+lpin=Bp: ， 转 步 又 b)。 上 文献 [17] 的 数值 实验 表明 使 用 前 者 的 效果 更 好 ， 因 此 ， 内 
定理 2 对 于 每 个 Ps ， wo 为 问题 (4) 的 一 个 稳定 点 。 令 “迭代 算法 停机 准则 设 定 为 

A 一 oo， 且 假设 严 为 fwe) 的 极限 点 ， 则 丈 为 均衡 约束 规划 问 上 YAOwe 咱 <E=10231 vr fw?) 

题 (3) 的 C- 稳 定点 。 其 中 : wo 为 外 迭代 算法 中 的 初始 点 。 同 时 ， 如 果 内 和 迭代 算 
证 明 Go) 入 (w) 的 定义 可 知 ，VG(W(VH() 为 第 ,ww 、 Vb a 

ij) 个 元 素 为 1， 其 他 位 置 元 素 为 0 的 向 量 。 则 法 先 代 一 次 结束 ， 则 减少 E 的 数值 ， 设 定 为 ==10=。 同 时 设 

{VG(W) (i e Zo (Ww)), VH ,(W) (j e Ty (w)), 定 p=1B=1.5 。 内 人 迭代 中 M=10。 设 定 最 大 外 迭 代步 数 为 50， 

VG,( 友 ),VH,(D (qe{q:G( 矶 >0,H,( 友 >0)} 线性 独立 。 由 文 ”最 大 内 迭代 步 数 为 200。 

献 [12] 的 引 理 3.2 得 ，w 为 互补 约束 规划 问题 (3) 的 可 行 点 。 本 文 首先 采用 100 个 随机 产生 的 绝对 值 方 程 问题 作为 算 

此 外 ， 又 因为 互补 约束 规划 线性 独立 约束 品 性 在 闷 成 立 ， 即 ”法 测试 问题 。 为 确保 问题 有 解 ， 即 使 得 4 的 奇异 值 大 于 1， 


{VG;(w),VHj(WD):ie To(w),je Tn( 罗 } 线性 独立 , 则 由 文献 [12] 的 定 本 文 取 4=Vzy ， 其 中 U 和 VV 为 均匀 分 布 在 [-10,10] 的 正 交 和 矩 


理 2.1 可 知 ， 丈 为 均衡 约束 规划 问题 (3) 的 C- 稳 定点 。 阵 ， 对 角 阵 的 元 素 从 区 间 [1,50] 中 随机 选取 。 取 均 匀 分 布 

在 外 迭代 算法 中 ， 本 文 每 进行 一 次 旬 代 都 需要 求解 步骤 ”在 [-1,1] 的 向 量 x* ， 则 5=4x-|x|。 外 迄 代 算 法 的 初始 值 

c) 中 的 问题 (4)。 下 面 给 出 内 迭代 算法 过 程 。 令 wo =|rand(2n,1)-rand(2n,D)| ， 本 文 算法 在 达到 设 定 的 误差 情形 下 
成 功 求解 了 所 有 测试 题 ， 具 体 实验 结果 如 表 1 所 示 。 


站 
0 表 1 100 个 随机 测试 是 计算 结果 


及 由 Gtw) 和 Hw) 的 定义 ， 罚 函数 问题 (4) 可 简化 为 如 下 表达 Table 1 Results of 100 random test problems 
形式 : 测试 题 结果 
min f(w) (5) 测试 题 维 数 n 1000 
st. w2>0. 算法 成 功 求解 个 数 100 
在 内 迭代 算法 中 ,本文 用 投影 Barzilai-Borwein 梯度 法 09 14x-1zxl- 吉 平均 值 3.6282e-007 
求解 问题 (5)。 搜 索 方向 设 定 为 外 从 代 步 数 平均 值 32.04 
d® =P(w® — ApaVf (WE) —w® x 一 x 平均 值 3.5772e-007 
其 中 : PO) 为 投影 算 子 ,Qs 为 由 BB 方法 获得 的 步 长 。 此 外 ， 同时 ， 本 文 也 采用 了 如 下 两 个 算 例 进行 数值 验证 。 
为 确保 全 局 收敛 性 ， 算 法 沿 上 述 投 影 方向 还 采用 了 非 单调 线 网 208] 和 矩 阵 4 定 义 为 
搜索 获得 线 搜索 步 长 4 ， 则 迭代 公式 为 weo =wo+Md0 。 di =4n,ain = a = =0.5,(i,j=1,2…n) ，D=(4-7De， 其 中 1 为 
为 区 别 内 外 迭代 算法 中 的 变量 记号 ， 本 文 令 内 迭代 的 变 mn 维 单位 矩阵 , e 为 元 素 全 为 1 的 列 向 量 。 该 问题 的 唯一 解 为 


量 。 潜 代 篇 尘 
:GD 。 初 始点 设 定 为 零 向 量 ， 算 法 对 =200.4=500, 及 
a) 设 置 初始 点 z? =we ，a 昌 = Vf(z5。。， 令 1=1; n=1000 的 维 数 问题 进行 测试 且 获 得 了 问题 的 唯一 解 ， 其 他 迭 
b) 如 果 1 PC -YfGcoD)-zol=0， 停 止 计算 ， 输 出 we =z2 ;人 代 次 数 , 停 机 时 的 罚 参 量 数值 及 14x-|x|-9l 数值 如 表 2 所 示 。 
否则 ， 令 4d0=PG0O-ogvyrdco)-z0，4=1, 转 步 又); 表 2 例 2 计算 结果 
c) 如 果 Table 2 Results of example 2 
ald ad) 维 数 n 人 迭代 步 数 罚 参量 P 上 4x-|x|- 可 
mag, fe) +107A( doy ve), 200 位 194.6195 2.7181E-007 
500 29 2.2168E+003 1.8910E-007 


则 令 zw =z?+H40 ， 转 步骤 d); 否则 令 2=4/4， 重 复 步骤 1000 29 3.3253E+003 6.2445E-007 
oj; 例 309] 和 抢 阵 4 定义 为 几 =500w =oi=1+randGtx 亡 。 其 中 


Im 
出 


录用 定稿 严 涛 : 
rand 为 [0,1] 的 随机 数 。 令 b=(4-7)e ,其 中 7 为 n 维 单位 矩阵 ， 
e 为 元 素 全 为 1 的 列 向 量 。 该 问题 的 唯一 解 为 x=(11…1)”。 


本 文 算法 对 n=200,n=500, 及 n=1000 维 数 问题 进行 测试 。 初 始 
点 设 为 零 向 量 。 经 过 有 限 次 迭代 ， 算 法 获得 问题 的 唯一 解 ， 
和 迭代 次 数 、 停 机 时 的 罚 参 量 数值 及 14x-|x|-9l 数值 如 表 3 所 
示 。 


基于 互补 约束 规划 模型 的 投影 


表 3 例 3 计算 结果 
Table 3 Results of Example 3 


维 数 n 和 帮 代 步 数 罚 参 量 P 1Ax-|x| -al 
200 20 129.7463 5.0775E-007 
500 25 656.8408 1.8032E-007 
1000 23 194.6195 4.7039E-008 

4 ”结束 语 


本 文 给 出 了 求解 绝对 值 方程 问题 基于 均衡 约束 规划 的 投 
影 梯度 型 算法 。 数 值 实验 验证 了 算法 的 有 效 性 ， 也 表明 本 文 
算法 为 求解 绝对 值 方 程 问 题 的 一 个 可 行 方法 。 
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